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Abstract 

We give a new procedure for generalized factorization and construction of 
the complete solution of strictly hyperbolic linear partial differential equations 
or strictly hyperbolic systems of such equations in the plane. This procedure 
generalizes the classical theory of Laplace transformations of second-order 
equations in the plane. 

1 Introduction 

Factorization of linear ordinary differential operators (LODOs) is often used in mod- 
ern algorithms for solution of the corresponding differential equations. In the last 20 
years numerous modifications and generalizations of algorithms for factorization of 
LODOs with rational function coefficients were given (see e.g J5])- Such algorithms 
have close relations with algorithms for computation of differential Galois groups 
and closed-form (Liouvillian) solutions of linear ordinary differential equations and 
systems of such equations (^H])- We have a nice and relatively simple theory of 
factorization of LODOs. 
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For linear j)arim/ differential operators (LPDOs) and the corresponding equations 
(LPDEs) the theory of factorization is much more difficult. To the best of our 
knowledge there are only a few theoretical results and only one algorithm for "naive" 
factorization of hyperbolic LPDO. In this introduction we will give a brief account 
of the previously obtained results and state our main result: existence of a recurrent 
procedure for non-trivial factorization and finding closed-form complete solutions of 
strictly hyperbolic systems of LPDEs in two independent variables with coefficients 
in an arbitrary differential field. 

Theoretically one may propose several very different definitions of factorization 
for LPDOs. The obvious "naive" definition suggests to represent a given opera- 
tor L = E 

|?|<m '^n - «n (•^)-^a;\-^x2 ■ ■ ■ -^x^ as a compositiou of lower-order LODOs: 
L = Li . . . Lk with coefficients in some fixed differential field. Unfortunately this 
definition does not enjoy good theoretical properties: a given LPDO may have sev- 
eral very different decompositions of this form, even the number of irreducible factors 
Ls may be different, as the following example (attributed in ^ to E. Landau) shows: 
if 

P = D^ + xDy, Q = D^ + 1, 
R = Dl + xD^Dy + ^, + (2 + x)Dy, ^ ' 

then L = QQp = RQ. On the other hand the second-order operator R is absolutely 
irreducible, i.e. one can not factor it into product of first-order operators with coef- 
ficients in any extension of Q(x, y). Still the "naive" definition of factorization may 
help to solve the corresponding LPDE in some cases; recently (|13j) an algorithm 
for such factorization for the case of hyperbolic LPDOs of arbitrary order was given. 

In 123] the adequate theoretical definition of factorization and a factorization 
algorithm for the case of overdetermined systems with finite-dimensional solution 
space and rational function coefficients was given. 

For a single second-order LPDO in two independent variables 

L = D^Dy - a{x, y)b^ - b{x, y)Dy - c{x, y) (2) 

we have a very old and powerful theory of Laplace transformations (not to be mixed 
with Laplace transforms!). We expose this nice theory in Section |21 Roughly speak- 
ing, an operator (j2)) is Laplace-factorizable if after several applications of differential 
substitutions [Dx- or £)j^-transformations), which change the coefficients of Q in a 
simple way, one obtains a naively-factorable operator L(fc) = {Dy + b(_k)){Dx + a(^k)) or 
Z/(fc) = {Dx + a(^k)){Dy + &(fc)). This phenomenon of non-trivial Laplace-factorization 
explains the existence of Landau example (^. The definition of Laplace-factorizable 
operators turns out to be very fruitful in applications, it was extensively used in clas- 
sical differential geometry (see e.g. j8;) and actively studied in the last decade in 
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the framework of the theory of integrable nonhnear partial differential equations 
PI |2ni • This is one of the most powerful methods of integration (construction of 
the complete solution with the necessary number of functional parameters) of the 
corresponding second-order equations in the plane. 

To the best of our knowledge the only serious effort to generalize the classical 
theory of Laplace-factorization to operators of higher order (in two independent 
variables) was undertaken in with rather obscure exposition but deep insight 
and a few enlightening remarks. Our approach, exposed in Section El gives a new 
uniform and general treatment of this topic directly for n x n strictly hyperbolic 
systems in two independent variables. Several modern papers jHl IE] investigate 
the theory of multidimensional conjugate nets initiated in jHl t. 4]; this line of 
research is in fact still in the domain of second-order equations in two independent 
variables: the systems discussed in the cited references are overdetermined systems 
with operators Q and solution spaces parameterized by functions of one variable. 
An interesting special case (operators (0) with matrix coefficients) was studied in 
|21L 122] , unfortunately the results are limited to this particular case of higher-order 
systems. 

A proper theoretical treatment of the factorization problem might be expected 
in the framework of the P-module theory (see e.g |^ and a very good exposition of 
the appropriate basic results in jTTj). Unfortunately even in this modern algebraic 
approach a "good" definition of factorization of LPDOs with properties similar to 
the properties of factorization of LODOs or commutative polynomials (decomposi- 
tion of algebraic varieties into irreducible components or primary decompositions in 
Noetherian commutative rings) is not an easy task. Without going into fine theo- 
retical details we refer to [22] where a variant of such "theoretically good" definition 
of generalized factorization of a single LPDO was given. As we have shown in [22], 
this definition generalizes the classical theory of Laplace-factorizable second-order 
operators. A drawback of this theoretical approach was lack of any factorization 
algorithm for a given LPDO. 

In the present paper we give a new procedure (generalized Laplace transforma- 
tions) for generalized factorization and integration of strictly hyperbolic LPDOs of 
arbitrary order with two independent variables or systems of such LPDOs. Section [H] 
is devoted to the detailed exposition of this new procedure. 

In Section ^ we give an example of application of this procedure to a 3 x 3 
system and construct its complete solution using the results of Section |21 After this 
a general scheme of generalized factorization and integration of a strictly hyperbolic 
system in the plane is given. We conjecture that this new procedure provides an 
algorithm for generalized factorization and closed-form complete solution precisely 
in the sense of if we limit the complexity of the answer. 
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2 The classical heritage: Laplace transformations 

Here we briefly sketch this classical theory in a slightly different form suitable for 
our purpose. The exhaustive exposition may be found in jHlCniEI- An arbitrary 
strictly hyperbolic second-order equation with two independent variables Lu = 
and the operator 

2 ^ ^ 

L = Y,PiDiDl'' + ai{x,y)D^ + a2ix,y)Dy + c{x,y), (3) 

1=0 

Pi = Pi{x, y), may be rewritten in characteristic form 

{X1X2 + aiXi + 0:2X2 + 03)^ = , . 

{X2X1 + aiXi + 02X2 + a3)u = 0, 

where ctj = ai{x,y), the coefficients of the first-order characteristic operators X^ = 
mi{x,y)Dx + ni{x,y)Dy are found (up to a rescaling Xj — > ji{x,y)Xi) from the 
characteristic equation m^po ~ "'^i^iPi + n^P2 = for the principal symbol of ©• 
Since the operators Xi do not commute we have to take into consideration in 
and everywhere below the commutation law 

[Xi, X2] = X1X2 - X2X1 = P{x, y)Xi + Q{x, y)X2. (5) 

Using the Laplace invariants of the operator (jlj): 

h = Xi{ai) + 0:102 - 03, k = X2(o2) + O1O2 - 03, 

we represent the original operator L in partially factorized form 

L = (Xi + 02) (X2 + ai)-h = (X2 + ai)(Xi + 02) - k. (6) 

From this form we see that the equation Lu = is equivalent to any of the first-order 
systems 

, . jX2U = -aiU + v, ^/e^ jXiU = -a2U + w, 
[ Xiv = nu — a2V. [ X2W = ku — aiw. 

Proposition 1 Any strictly hyperbolic LPDE is equivalent to a 2 x 2 first-order 
characteristic system 

{ Xiui = aii{x,y)ui + ai2{x,y)u2, 

\ X2M2 = 02l(x,?/) Ml + 022(0;, 2/)u2, 

with Xi = mi{x,y)Dx + ni{x,y)Dy, Xi 7^ ^{x,y)X2, and any such system with 
non-diagonal matrix (o^) is equivalent to a second-order strictly hyperbolic LPDE. 
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Proof. Transformation of a strictly hyperbolic LPDE into the form (jSI) is already 
given. The converse transformation is also simple: if for example 0:12 7^ then 
substitute U2 = {XiUi — aiiUi)/ai2 into the second equation of the system (jSJ. □ 

Proposition 2 // a 2 x 2 first-order system 

vA^faii ai2\fvi\fbu 6i2Wfi\ 

V2j^ \a2l a22 ) V&2I ^22/ \V2 ) 

with aij = aij{x,y), h-ij = bij{x,y) is strictly hyperbolic (i.e. the eigenvalues Xk{x,y) 
of the matrix (aij) are real and distinct), then it may be transformed into a system 
in characteristic form (0) ■ 

Proof. Let Ai(x,|/), A2(a;,?/) be the eigenvalues of (ajj) and pi = {pii{x,y),pi2{x,y)), 

P2 = {p2i{x,y),P22{x,y)) be the corresponding left eigenvectors: J2k 

Form the operators Xj = — XiDy and the new characteristic functions 

Ui = HkPikVk- Then XiUi = J2k{^iPik)vk + HkPik {{vk)x - \{vk)y) = J2k,sPik{aks - 

KSks){vs)y+J2k,sPikbksVs+J2k{XiPik)vk=J2sVs {j2kPikbks+ {XiPis)'j = "Ek UkCXik{x , y) , 

so we obtain the characteristic system (jSI). □ 

The characteristic system (jH)), equivalent to 0, is determined uniquely up to 
operator rescaling Xi — > 'yi{x,y)Xi and gauge transformations Ui — > gi{x,y)ui. It is 
easy to check that the gauge transformations to not change the Laplace invariants 
h = X2(aii) - Xi (022) - X1X2 ln(Q;i2) - Xi{P) + Pan + ai2a2i + (^22 + X2 (In 012) + 
P)Q and k = ai2«2i5 they are just the Laplace invariants of the operator (0)), 
obtained after elimination of U2 from (jSJ. Rescaling transformations of Xj change 
them multiplicatively h 7172^, k 7172^- From the proofs we see that for a fixed 
equation Lu = with the operator (jS)) we obtain two different (inequivalent w.r.t. the 
scaling and gauge transformations) characteristic systems (|7j) and from every fixed 
system (jH)) we obtain two different (inequivalent w.r.t. the gauge transformation 
u —>■ g{x,y)u) hyperbolic LPDEs: one for the function Ui and the other for the 
function U2- This observation gives rise to the Laplace cascade method of integration 
of strictly hyperbolic LPDEs Lu = with operators Q: 

(£1) If at least one of the Laplace invariants h or k vanishes then the operator L 
factors (in the "naive" way) into composition of two first-order operators as we see 
from if we perform an appropriate change of coordinates (x, y) (x, y) (NOTE: 
for this we have to solve first-order nonlinear ODEs dy/dx = ni{x,y)/mi{x,y), 
cf. Appendix in jTH]) one can suppose Xi = D^, X2 = Dy so we obtain the 
complete solution of the original equation in quadratures: if for example Lu = 
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0^ + a2{x, y)){Dy + ai{x, y))u = 0, then 

u = exp ^— J ai dy^ (^X(x) + J Y{y) exp (^J (ai dy — dx)^ dy 
where X{x) and Y{y) are two arbitrary functions of the characteristic variables x, 

y- 

(£2) If /i 7^ 0, A; 7^ 0, transform the equation into one of the systems ^ (to fix 
the notations we choose the left system {Si}) and then finding 

u = (Xif + a2v)/h (10) 

substitute this expressions into the first equation of the left system (^i) in ((Tj), 
obtaining a Xi -trans/ormed equation = 0. It has Laplace invariants (cf. ^J) 

/i(i)=Xi(2ai-P)-X2(a2)-^i^2ln/i + QX2ln/i-a3+(ai-P)(a2-Q) 
= 2h-k-XiX2\nh + QX2\nh + X2iQ)-MP) + '^PQ, (11) 
= h. 

If = 0, we solve this new equation in quadratures and using the same differential 
substitution (fTUj) we obtain the complete solution of the original equation Lu = 0. 

(£3) If again 7^ 0, apply this Xi-transformation several times, obtaining a 
sequence of second-order operators 1/(2), 1/(3), ... of the form (jlj). If on any step we 
get /i(fc) = 0, we solve the corresponding equation L(^k)U(k) = in quadratures and, 
using the differential substitutions (fTUj). obtain the complete solution of the original 
equation. Alternatively one may perform X2-transformations: rewrite the original 
equation in the form of the right system (^2) in ((Tj) and using the substitution 
u = {^(2^ + aiw)/k obtain the equation Z/(_i)-u; = with Laplace invariants 

/i(-i) = fc, 

... (12) 
fc(_i) = 2k-h~ X2X1 Ink- PXi \nk + X2{Q) - Xi{P) + 2PQ. 

In fact this X2-transformation is a reverse of the Xi-transformation up to a gauge 
transformation (see y|^). So we have (infinite in general) chain of second-order 
operators 

X2 f X2 f X2 f ^1 f -^^ f -^^ f o\ 

. . . ^ -t^(-2) ^ -^(-1) ^ -t>(i) L(2) ^ . . . {16) 

and the corresponding chain of Laplace invariants 

...,/l(_3), /i(-2), ^(-1), ho = h, /l(i), /l(2), /l(3),--- (14) 
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with recurrence formulas (jllj) . (jl2j) . We do not need to keep the invariants in 
(HD) since A;(j) = If on any step we have /i(Ar) = then the chains (fT^ and 

(HH) can not be continued: the differential substitution (fTUI) is not defined; precisely 
on this step the corresponding LPDE is trivially factorable and we can find the 
complete solution for any of the operators of the chain (fT^ . For simplicity let us 
choose characteristic variables {x,y), so that Xi = D^, X2 = Dy. The complete 
solution of the original equation in this case has the form 

u = co(x,y) (F+JGpdy) + c^(x,y) (F'+JF§dy) + 

... + c„{x,y){F(^^ + JGffdy). 

where F{x), G{y) are two arbitrary functions of the characteristic variables and 
Ci(x,y), P(x,y) are some definite functions obtained in the process of Laplace trans- 
formations from the coefficients of the operator Q. As one may prove (see e.g. |H]) 
if the chain (fT^ is finite in both directions (i.e. we have /i(Ar) = 0, /^(-x) = for 
some > 0, > 0) one may obtain a quadrature-free expression of the general 
solution of the original equation: 

u =coF + ciF' +. . . + cjvF(^)+ doG + diG' + . . . -I- dK+i&^^^^ (16) 

with definite Ci{x,y), di{x,y) and F{x), G{y) — two arbitrary functions of the 
characteristic variables and vice versa: existence of (a priori not complete) solution 
of the form (|16|) with arbitrary functions F, G of characteristic variables implies 
h(^s) = 0, = for some s < N, r < K. So minimal differential complexity of 

the answer ()16|) (number of terms in it) is equal to the number of steps necessary 
to obtain vanishing Laplace invariants in the chains p3|) . (fT^ and consequently 
naively-factorable operators. Complete proofs of these statement may be found in 
P t. 2], UnmH for the case Xi = D^, X2 = Dy, for the general case cf. [m p. 30] 
and n. 

Example 1. As a straightforward computation shows, for the equation u^y — 
^{x+i/f M = the chain ()14|) is symmetric {h(i) = /i(-i-i)) and has length n in either 
direction. So the complexity of the answer ()16|) may be very high and depends on 
some arithmetic properties of the coefficients of the operator Q; for the equation 



Uxy— (x+y)'^ '^ ~ ^ chains (fT^ . (fT^ will be infinite unless the constant c = n(n+l). 

Example 2. For a stochastic ODE x = p{x) + a{t)q{x) with binary (dichotomic) 
noise a(t) = ±1 and switching frequency u > the averages W{x,t) = {W{x,t)) 
and Wi{x,t) = {a(t)W{x,t)) for the probability density W{x,t) in the space of 
possible trajectories x{t) of the ODE satisfy a system of the form Q (see [T^): 

Wt + (p(x)iy), + iqix)W,), = 0, 

{W,)t + 2uWi + {p{x)W,)^ + {q{x)W)^ = 0. 
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The characteristic operators and left eigenvectors are simple: Xi = Df — XiD^, 
Ai^2 = —p{x) ± q{x), pu = P21 = P22 = 1, P12 = — 1- The characteristic system (jS)) 
for the new characteristic functions ui = W — Wi, U2 = W + Wi is 

XiUi = -{p^ - + v) ui + V U2, ^^g^| 
X2U2 = I' ui - + + u) U2. 

The Laplace invariants are h = v"^ — [pxxq^ip + q) +Pxl'^ ~ Pxqxqi^^P + q) — qxxPqip + 
q)qlp{2p + q)]/q'^, k = z/^, so if u, p{x) and q{x) satisfy a second-order differential 
relation h = 0, one can solve (fTTjl in quadratures. Especially simple formulas may 
be obtained for polynomial p{x) = pix + P2X^, q{x) = q2X^, pi > 0, p2 < 0: in this 
case k = u"^, h = /i(-2) = v'^ — p\ so ii v = pi after the necessary transformations 
we obtain the following quadrature-free expression for the complete solution of the 
system (fT7|) : 



xG'{y) + pI{p2X + pi)G{y)\ , 

where x = —t + — -r-^- — ^, y = —t + — In — -j^ — ^ are the characteristic 

variables {X2X = 0, Xiy = 0) and F, G are two arbitrary functions of the cor- 
responding characteristic variables. For the case z/^ 7^ pf we can compute other 
Laplace invariants of the chain (|T^ : h(i) = h{^-z) = — 4^^, h(2) = /i{-4) = u"^ — 9p1, 
h(3) = ^(-5) = — 16j9^, ... so for the fixed p{x) = pix + p2x'^, q{x) = q2X^ and 
z/ = ±pi, V = ±2pi, V = ±3pi, ... one can obtain closed-form quadrature- free 
complete solution of the system (fT7j) . with increasing complexity of the answer (fTBj) . 

Remark. The forms (|15j) . for the complete solution are local: in the general 
case due to nontrivial topological picture of the trajectories of the vector fields Xi, 
X2 in the plane we are unable to guarantee existence of the global coordinate change 
{x,y) I—*- {x,y) (cf. Example 2 above). 

3 Generalized Laplace transformations of n x n 
hyperbolic systems 

Hereafter we suppose that the LPDE of order n > 2 

Lu= PiA^,y)D'xDiu = (19) 

i+j<n 



w 



Q2 



x^ 



F'{x)-p,Fix)+p,G'{y) -piG{y) 



-q2xF'{x) + pi{p2X + pi)F{x) + piq 
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is strictly hyperbolic, i.e. the characteristic equation ^ Pi,j^^ = has n simple real 

i+j=n 

roots Xk{x,y). 

Proposition 3 Any strictly hyperbolic LPDE ^19^ is equivalent to a nxn first-order 
system in characteristic form 

XiUi = aik{x, y)uk. (20) 

k 

Proof. The principal (rath-order) part of (|T9|l decomposes into the product of the 
characteristic operators Xi = — XiDy modulo lower-order terms. Other lower- 
order terms of any given order s may be also written as sums of products of Xi 
(modulo terms of order < s) in the following unique way: 

n-l 

L = XiX2---Xn + Yl (^s,h-is{x,y)Xij^- ■ ■Xi^+ ao{x,y). (21) 

s=l ii<...<is<s+l 

This is easily proved using induction over s. For s = 1 the terms of order 1 are 
Pifi{x,y)Dx + Po^i{x,y)Dy = ai{x,y)Xi + a2{x,y)X2 since the eigenvalues Ai, A2 are 
distinct. If this is proved for terms of order < si we may suppose (after a local 
change of variables) that Xsj^+2 = so ^ pijDlD^y = po,si+i-D*^^^+ 

i+j=si+l 



Y: P^+l,DlDl = qix, y)X,X2 ■ ■ ■ + ^ P.+,,DIDI = 

yi+j=si j \i+i=sx ) 

q{x, y)XiX2 ■ ■ ■ ^si+i + Yl (^suh-is^ i^, y)^h ■ ■ ■ ^is^ (modulo terms 

\ii<...<jsi<si+l / 

of lower order) . 

Now, provided we have the form (PT|) of the operator L, introduce the charac- 
teristic functions m„ = m, Un-i = X„m„ + Un-2 = Xn-iu^-i + a„_i,„-iM„_i + 
= X2U2 + YJl=2 o.2sUs with SO far undefined y). We 
have to define them inductively so that the expression M = XiUi + Z]s=i ctis^s 
(again with ais{x, y) to be determined) would coincide with the operator L. In fact, 
the highest-order term in M is just X1X2 ■ ■ ■ The terms of order n — 1 in M are 
just n products Xj^ ■ ■ ■ ii < . . . < in-i < n (in every such product exactly one 
of Xj, 1 < i < 72, is missing) with coefficients atk-, 1 < A; < n, so we just fix atk to 
be equal the respective an~i/i^-i^_^ in (0T|) . This makes M — L to be of order n — 2. 
Each coefficient at the terms Xj^ ■ ■ ■ Xj„_2 of order ra — 2 in M — L includes deriva- 
tives of the already defined akk and precisely one of ak-i,k, 2 < k < n, so choosing 
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the latter appropriately we lower the order of M — L again. This step-down process 
will give us all apg so we obtain the following system, equivalent to Lu = 0: 

X2U2 = ~YTs=2(^2sUs + Ui, 

□ 

Remark 1. There are many ways of transforming a given higher-order LPDE 
into a first-order system. But for our purposes we need a precise equivalence, guar- 
anteeing a one-to-one correspondence between the solution space of the original 
equation and the solution space of the system. If for example we take the equation 
Uxx — Uyy — a{x, y)u = and (after [71 Ch. I]) form the following first-order system 

Ux = P, 

' Qx = Py, (22) 
, Px = qy + a{x,y)u, 

then we shall fix for it limited Cauchy data q\x=o = Uy\x=o (so taking some subset of 
solutions of (1221)) order to have a one-to-one correspondence between the solution 
space of the original equation and some subset of the solution space of ((221) • Another 
straightforward possibility 

' Ux=p, 

< Uy = q, (23) 
^ Px = qy + a{x,y)u 

gives in fact an overdetermined degenerate system, it can not be put into the stan- 
dard form (^^ below. 

Remark 2. The converse of Proposition El is not true: already for = 3 a generic 
system (l2Uj) gives for each Ui an overdetermined involutive system of equations (see 
for example T, Ch. I]). So hereafter we will construct our generalized Laplace trans- 
formations directly for systems in characteristic form (j^Uj) with Xi ^ fi{x,y)Xj for 
i ^ j- But for n X n first-order systems the equivalence can be proved. This is a 
well-known fact (see for example [Zj), we give its simple proof for completeness of 
our exposition. 

Proposition 4 Any n x n first-order linear system 

n n 

{vi)x = (^ik{x, y) {vk)y + 51 ^ik{x, y)vk (24) 

k=l k=l 
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with strictly hyperbolic matrix {aik) (i.e. with real and distinct eigenvalues of this 
matrix) is equivalent to a system in characteristic form (ji^^) . 

Proof. A straightforward calculation similar to that in the proof of 
Proposition 121 shows that if pi = 

{pii{x,y), . . .,pin{x,y)) are the left eigenvectors, 'EkPikakj = Kpij, Xi = - XiDy 
are the characteristic operators, then the new characteristic functions Ui = J2k Pik^k 
satisfy a system of the form (j2Uj) . □ 

The characteristic form ()20j) of a given system ()24|) is defined up to operator 
rescalings Xj — 'yi{x,y)Xi and gauge transformations Ui gi{x,y)ui. The quanti- 
ties = 

^ 7^ J 7^ ^5 ^iicl other similar cyclic products of 
the coefficients of (|2Up are invariant w.r.t. the gauge transformations so it is natural 
to call them Laplace invariants of ()20p . The complete set of such invariants can be 
found using the standard methods [iSilll- We will not dwell upon this problem. 

Our generalization of the Laplace transformations is based on the following ob- 
servation. For a given 2x2 system (jH)) the procedure of Xi-transformation consists 
of the following steps: 

* / X2U = -aiu + v f 

[ XiV = hu — Oi2V 

If we will try to reformulate this in terms of 2 x 2 characteristic systems (jH)), we 
arrive at the following procedure: 

(£1) For a given system (jHJ choose one of the two variables Ui and eliminate it, for 
example we may eliminate U2 using the first equation of (jHl): M2 = (^1^1— oii^i) / ai2- 
Substitute it into the other equation of the system obtaining a second-order equation 
for uy. X2(^(XiUi - anWi)) - asi'Ui - ^(XiUi - anWi) = 0. 

(£2) Swap the operators X\ and X2 in the second-order term of this equa- 
tion using the appropriate commutation laws and then rearrange the first-order 
terms obtaining an equal expression of the form Xx {-^[X2U\)^ -|- pi(a;, -(- 

P2[x,y)X2Ux^P'>,[x,y)ux = Xi (^(XgWi + piai2Ui)^+^{X2Ui+piai2Ui)+p^{x,y)ui = 
0. 

(£3) The last form of this second-order equation immediately suggests to set U2 = 
{X2U1 + piOuUi) / ai2 so we can rewrite the equation in the form of the transformed 
2x2 characteristic system 

X2U1 = -pi«i2Mi + ai2M2, ^25) 
X1U2 = -P2U1 -p^{x,y)u2. 

Remark. One should not be mistaken: if we apply this "swapping procedure" twice 
we do not necessarily return to the original system. This depends on the choice 
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of the characteristic function to ehminate on step (£i) in (|25|): ehminating U2 we 
certainly return to the original system (jS)), but eliminating Ui from (|25p we actually 
move one more step forward in the chain (fT^ . 

This commutation trick gives the basis for our generalized Laplace transforma- 
tions: 

(Ci) For a given nxn characteristic system (j2(J|) choose one of its equations with 
a non- vanishing off-diagonal coefficient a^fc 7^ 0, find Uk = {XiUi — J2sj^k(^isUs)/c(ik 
and substitute this expression into all other equations of the system. We obtain one 
second-order equation 

Xk{ — {XiUi - ^^aisUs)) -Yl (^kpUp - —{XiUi -Y aisUs)=0 (26) 

•^^fc sjtk p=ik 

and n — 2 first-order equations 

^3^3 ~ Y '^jsUs - —{XiUi - Y o^isUs) = 0, j ^ i, k. (27) 

S^fc g-^fc 

(£2) The second step consists in rewriting the system (j2SI); in the following form 
with slightly modified unknown functions Uj = Uj + Pj{x, y)ui, j 7^ i, k, Ui = Ui, new 
coefficients Ppg{x,y) but the same characteristic operators Xp-. 

Xi{ — {XkUi - JZ A^^J) - Yl^kpUp - —{XkUi A^Ws) =0, (28) 

s^k p^k s^k 

^j'^j - ~ —{XkUi - Y l^isUs) =0, i ^ i, k. (29) 

sj^k 

As we prove below this is always possible in a unique way. 

(£3) Introducing Uk = -^{XkUi — J2s=ik PisUs) rewrite (j28|) . (j29|) as the trans- 
formed characteristic system 

{X-lUk ^2p f-^kpUipj 
XkUi = J2sjtk PisUs + CXikUk, (30) 
XjUj = EsPjsUs, i^i.k. 

Correctness of the step {€2). First of all we rewrite the l.h.s. of ()26|) (again using 
the appropriate commutation laws) as Ei = Xj(— (X^-Uj)) + piXiUi + pkXkUi — 

'Esyti,k^^kUs + Hs^kCisUs and the l.h.s. of (j2II) as Ej = XjUj - Y^s^kCjsUs - 
PjkXkUi — pjjXjUi, where we have used the fact that any 3 characteristic operators 
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are linearly dependent (we have only 2 independent variables!): Xj = (j)j{x,y)Xk + 
ipj{x,y)Xj. Introducing new Uj = Uj — pjj{x,y)ui, j ^ i,k, Ui = Ui we have 
Ej = XjUj - J2s^kCjsUs - PjkXkUi and Ei = Xi{-^{XkUi)) + piXiUi +p^,XkUi - 

J2s^i,k ^XkUs + Es^k CisUs = Xi{^{XkUi)) +PiX,Ui+PkXkUi - Y.s^i,k J^i^iUs - 

ipsXsUs) + Hs^kCisUs- So adding to it a suitable combination of Ej we obtain 
E, = Ei- Jff^^s = Xi{^^{XkUi + a,kPiU^ - ^u,)) + (^(X^m, + 

c^ikPiUi — Y^sj^i.k ^"^s)) + Y.sjtk'^isUs, which coiucides with the basic coefficients 
Pks{x, y), Pi{x, y) are therefore defined. Now is in fact the equation Ej = with 
slightly rearranged terms. □ 



4 An example and the general factorization pro- 
cedure 

Let us consider the following 3x3 system 

XiUi = Ui + 2U2 + Ms, 

< X2U2 = -Qui + U2 + 2^3, (31) 

X3U3 = 12Mi + 6U2 + Ms, 

with Xi = Dx, X2 = Dy, X3 = Dx + Dy. Let us eliminate M3 using the ffist equation: 

M3 = XiUi — Ml — 2u2. (32) 

Substituting this expression into the other equations we obtain the following two 
equations: 

E2 = X2U2 + Qui — U2 — 2(Xini — Ml — 2M2) = 0, 

E3 = X3(XiMi — Ml — 2^2) — 12Mi — 6M2 — (^iMi — Ml — 2M2) = 0. 

This completes the step {Ci) of our generalized Laplace transformation. 
Step (£2) is done in 4 substeps: 

1) Change XiMi to X3M1 — X2M1 in E2 and collect the terms with X2: E2 = 
X2U2 + 6mi — M2 — 2X1M1 + 2mi + 4m2 = X2M2 + 8mi + 3m2 — 2X3M1 + 2X2M1 = 
-^2(m2 + 2mi) + 8mi + 3m2 — 2X3M1 = 0, so we now can introduce the new function 
M2 = M2 + 2mi. Change everywhere M2 to M2 ~ 2mi: 

E2 = X2U2 + 2mi + 3m2 - 2X3M1 = 0, 

E3 = X3X1M1 + 3X3M1 - 2X3M2 - 3mi - 4m2 - XiMi = 0. 

2) Swap X3 and Xi in E^, substitute X3M2 = X1M2 + X2M2 and collect the terms 
with Xi. 
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i?3 = Xi{X3Ui — Ml — 2U2) + SX3U1 — 2X2^2 — 3mi — 4^2 = 0. This suggests to set 
U3 = {X3U1 -ui- 2U2). 

3) Change everywhere X^ui to + Ui + 2^2: 

E2 = X2U2 + 2Mi + 3U2 - 2(m3 + Ml + 2112) = X2li2 - U2 - 2II3 = 0, 

= X1U3 + 3{u3 + ui + 2U2) - 2X2U2 - 3ni - 4^2 = X1M3 + 3u3 - 2X2U2 + 2^2 = 0. 

4) Now get rid of the term —2X2^2 in E^: 
E^ = E; + 2E2 = X1U3 -^3 = 0. 

The resulting equations E2 = X2U2 — U2 — 2m3 = 0, E3 = X1U3 — U3 = and 
U3 = (X3M1 — Ml — 2M2) give us the transformed system: 

X1II3 = M3, 

< X2M2 = 2^3 + U2, (33) 

X3M1 =U3 + 2U2 + Ml. 

Its equations can be consecutively solved: the first equations gives II3 = e^F{y). 
Using the standard method of the variation of constants we solve the second equa- 
tion obtaining U2 = e^(G'(x) + 2/ e^~^F{y) dy). In order to remove the quadrature 
we introduce a new arbitrary function F{y) instead of F{y): F{y) = e^F'{y). In- 
tegrating by parts we obtain U2 = e^{G{x) + 2e^F{y)). The last equation gives 
Ml = C(x, |/)exp^^ where C{x,y) satisfies the equation X3{C) = exp ~^^^ {u3 + 
2U2) = exp ^^{F'{y) + AF{y)) + 2exp^^G{x). Again introduce new arbitrary 
functions G{x), F(y), such that G(x) = G'{x), F{y) = e~'^F'{y), then we can find 
G{x, y) = 2 exp ^G{x) + exp ^{3F{y) + T {y)) + H{x - y). Finally we get the 
following complete solution of the transformed system (jHHj) : 

' Ml = 2e^G(x) + e^(3F(?/) + F\y)) + exp ^H{x - y), 
< M2 = eS'G'(x) +2e^'F'(?/), 
^ U3 = e%F"{y)-F'{y)), 

where F{y), G{x) and H{x — y) are three arbitrary functions of one variable each. 
The complete solution of the original system (p?T|) can be easily found using the 
differential substitutions of the steps (£1) and {C2): M2 = M2 — 2mi, M3 = XiMi — 
Ml — 2m2 = XiMi + 3mi — 2m2 (and certainly mi = mi!). 
The general factorization procedure. 

I) If we have to solve an equation Lm = or a system (j24|l transform it into the 
characteristic form ()2()|1 . 

II) If the matrix {aij{x,y)) of the characteristic system is upper- or lower- 
triangular (similar to the matrix of ()33p) solve the equations consecutively. 



14 



III) If the matrix is block-triangular, the system factors into several lower-order 
systems; try for each subsystem the step IV. 

IV) In the general case of not-triangular matrix {aij{x,y)) perform several (con- 
secutive) generalized Laplace transformations, using different choices of the pivot 
element aik 7^ on step (£1). The goal is to obtain a block-triangular matrix for 
one of the transformed systems. 

The main problem of this procedure is absence of upper bounds for the number 
of steps. As we have seen in Section |21 {Example 1) this bound depends on the 
coefficients of the equation (system) in a nontrivial way. On the other hand in this 
situation (the classical Laplace-factorization and our generalization) the procedure 
of transformation does not require solution of differential equations and does not 
depend on the differential field of coefficients, contrary to the case of LODOs. 
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